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A neutral quantum particle with magnetic moment encircling a static electric charge acquires a 
quantum mechanical phase (Aharonov-Casher effect). In superconducting electronics the neutral 
particle becomes a fluxon that moves around superconducting islands connected by Josephson junc- 
tions. The full understanding of this effect in systems of many junctions is crucial for the design 
of novel quantum circuits. Here we present measurements and quantitative analysis of fluxon in- 
terference patterns in a six Josephson junction chain. In this multi-junction circuit the fluxon can 
encircle any combination of charges on five superconducting islands, resulting in a complex pattern. 
We compare the experimental results with predictions of a simplified model that treats fluxons as in- 
dependent excitations and with the results of the full diagonalization of the quantum problem. Our 
results demonstrate the accuracy of the fluxon interference description and the quantum coherence 
of these arrays. 



The formation of macroscopically large coherent 
states in systems with many unquenched degrees of 
freedom tests our understanding of quantum mechan- 
ics and it is essential for quantum computation. One 
of the most striking consequences of such coherence 
are interference patterns that are expected to ap- 
pear when a charged particle encircles a magnetic flux 
(Aharonov-Bohm effect [1J) or when a flux encircles 
a charge (Aharonov-Casher effect [2]). The quantum 
coherence implied by these effects is a fragile phe- 
nomenon which is easily destroyed by uncontrolled de- 
grees of freedom. In an ideal Josephson junction array 
most microscopic degrees of freedom are quenched by 
electron pairing into Cooper pairs; the only remaining 
degrees of freedom are the phase of the order parame- 
ter on each island or the charge conjugated to it. Ob- 
servation of the interference provides the evidence of 
the full control of the quantum system in the ground 
state. In Josephson junction arrays it proves the ir- 
relevance of the uncontrolled degrees of freedom such 
as two level systems, non-equilibrium quasi-particles, 
etc. 

Aharonov-Bohm (AB) effect in small and large 
Josephson arrays is a very well established phe- 
nomenon: in the former it leads to critical current 
oscillations in SQUIDs|3j, in the latter it results in 
a complicated magnetic field dependence with many 
peaks at commensurate fields [IH7]- The experimental 
confirmation of its dual, the Aharonov-Casher (AC) 
effect, is less clear. It was observed for small Joseph- 
son circuits where vortices moved in a ring encircling 
a single charge |Sj. However, large arrays studied in a 
number of works show the appearance of the interme- 
diate "normal" phase of the arrays which is character- 
ized by a non-zero resistance |SJ O [TO]. Non-zero re- 
sistance implies that the fluxon motion is dissipative; 
this excludes quantum coherence. It is very important 
to establish the presence or absence of this dissipation 
and its possible origin in well controlled medium size 



arrays. This is the main goal of our work. 

The duality of AB and AC effects can be illustrated 
by analyzing the quantum mechanical phase resulting 
from the braiding of particle with charge q and a neu- 
tral particle with magnetic moment /2. If the charged 
particle is at rest while the neutral one moves, the 
former generates an electrical field at the posi- 
tion of the latter that gives the interaction energy 
1 1 = x E^] ■ Conversely, the magnetic mo- 
ment at rest generates a vector potential A q at the 
position R q of the moving charge, that gives the in- 
teraction energy I2 = - A q ■ R q . In either case, the 
acquired phase is given by the time integral of the 
interaction energy. In case of the moving charge, 
this phase is 8<\>ab = \jh)§Aq ■ dR q (AB effect); 
in case of a moving magnetic moment this phase is 

5((>ac = 7^ § [P x E^j ■ <iR M (AC effect). Experimen- 
tally the former was first observed 50 years ago as 
an electron interference pattern in magnetic field [ 1 1 J : 
the latter was measured with percent accuracy in neu- 
tron and atomic interferometry experiments |12[ 113] . 
In Lorentz invariant systems of neutral and charged 
particles the distinction between the two effects is im- 
possible. What seems as AB effect for the observer 
in the rest frame of the neutral particle becomes AC 
effect for the observer in the rest frame of the charged 
particle. In solid state devices these effects are distin- 
guishable because the rest frame is fixed by the device, 
and therefore the observation of AB interference does 
not imply the one of AC and vice versa. Because un- 
controlled degrees of freedom turn out to be mostly 
electric charges (either background charges or non- 
equilibrium quasi-particles in superconductors), ex- 
perimentally it is difficult to remain in the rest frame 
of the charge and, consequently, the observation of the 
AC effect is much more challenging. 

A Josephson junction circuit can be described ei- 



ther in terms of the superconducting phase or in terms 
of the charges of its islands. If the charging energy 
Ec — e 2 /(2C) is larger than the Josephson energy 
Ej, Cooper pairs arc almost localized. The dynam- 
ics of a Josephson junction array can be viewed as 
due to the rare motion of these pairs. In the oppo- 
site limit Ej ^> Ec the charge is delocalized and the 
array dynamics can be viewed as due to rare phase 
slips resulting in the motion of fluxons. The nature of 
the elementary excitations does not preclude however 
the description of the circuit in the charge (or phase) 
basis. The elementary excitations only become more 
complicated objects for Ej > Ec in the charge basis 
or for Ec > Ej in the phase basis. 

In the following, we analyze the ground state prop- 
erties of a 6 Josephson junction chain as a function 
of the gate voltage that induces polarization charges 
for Ej/Ec = 2 — 3. Because in this regime the in- 
dividual excitations are fluxons, the properties of the 
circuit are due to the interference between phase-slips 
on different junctions. The interference pattern is due 
to the charges induced on the array islands which is 
exactly AC effect. The difference between the longer 
chain (of six junctions) studied here and the previ- 
ous works[Hl ES] is that fluxons can take one of the 
six possible routes resulting in a much more compli- 
cated interference. In the following we compute the 
expected properties of the array assuming that phase 
slips are independent excitations. Because this as- 
sumption might be questioned for Ej/Eq = 2 — 3 
we have also performed the complete diagonalization 
of the Hamiltonian. Finally we compare the results 
of both approaches to the measured data. Our main 
conclusion is that the phase slip approximation pro- 
vides a semi-quantitative description of the data in 
this regime and that the observed interference pat- 
tern evidences the quantum coherent dynamics of our 
relatively large circuits. 

Figjlji shows an idealized view of our circuit: a su- 
perconducting ring containing five islands connected 
by Josephson junctions. A gate voltage Vc induces 
the charge frustration qi — CfVc/(2e) + q® on the 
i-th island. Here q® are static offset charges. The 
couplings to the gate electrode Cf are not equal for 
all islands and induce a general charge configuration 
K — Q2i 93j 94; 95), expressed in units of 2e. White 
traces in Fig[T^ represent six possible paths for a 
fluxon to cross the ring, through one of the six Joseph- 
son junctions. The probability of this event is given by 
the quantum phase-slip amplitude vj of a single junc- 
tion. Vj contains an AC phase-factor depending on the 
islands charges qi (see eq. ([TJ and Q). The ground 
state of the SQUID chain depends on the Coherent 
Quantum Phase-Slip (CQPS) amplitude that results 
from the macroscopic interference of six fluxons. The 
CQPS amplitude, v* , is obtained by summing up all 
phase-slip amplitudes on the junctions [15 (for com- 
putation details see the Supplementary Information 
text): 




Figure 1: Schematic view of the experimental setup used to 
probe the phase-slip interference in a chain of 6 Josephson 
junctions. In (a) we show an idealized view of the exper- 
imental design. The chain contains five small supercon- 
ducting islands connected to each other and to the leads 
by identical Josephson junctions. The islands are coupled 
to a nearby gate electrode. In (b) we present the electri- 
cal scheme of the measurement. The six-SQUID chain is 
inserted in a superconducting loop. The flux $c created 
by on-chip coils controls the phase difference 7 over the 
chain. The independently controlled flux $s through the 
SQUID loops is used to tune in situ the Josephson cou- 
pling Ej = E°j cos (7r<&s/<i>o), where E° = 2K and $0 is 
the magnetic flux quantum. The charging of one SQUID 
is Ec = 660 mK. The phase difference over the read-out 
junction is denoted by S. The gate electrode couples to the 
charge qi on island i via the capacitance Cf. The coupling 
to the central island Cf is at least 10 times larger than 
all other capacitances and determines the dominant gate 
effect at low voltage, (c) shows the calculated ground and 
first excited state for the 6 junction chain as a function of 
the phase bias 7 and the induced charge 53 on the central 
island, for charge configurations of the type (0, 0, 53, 0, 0). 



Here v is the magnitude of the phase-slip amplitude 
for a single Josephson junction. In the quasi-classical 
approximation, valid at Ej 3> E c , it is|16|: 
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The first two energy levels of a Josephson junction 
chain are shown in Fig. [IJ. These energy levels have 
been calculated by diagonalizing the Matveev-Larkin- 
Glazman (MLG) tight-binding Hamiltonian [15J: 



H \m) = E rn \m) - v* [|m- 1) + \m + 1)] (3) 

Here E m = ~ 2irm) 2 is the energy of the |m) 

state of the chain polarized at phase 7 and m is the 
quantum variable that counts the number of vortices 
having crossed the chain through one of the junctions. 
The model ([3| makes two important assumptions: the 
quantum phase slips on different junctions lead to the 
same quantum states and these tunneling processes 
are independent events. As it can be seen from eq. 
Q, the AC interference of CQPS is an intrinsically 
2e periodic effect. 

In our sample, each junction of the chain is real- 
ized by a SQUID (see Fig. [TJd) to enable tunable 
Josephson coupling Ej. Consequently we can control 
in-situ the strength of the quantum phase slip am- 
plitude v through the magnetic flux $5. To measure 
the CQPS effect on the ground state of a Josephson 
junction chain, we have shunted the chain by a large 
read-out junction (see Fig. [Tja, [T7j and (TH] ) . The 
flux $c m the superconducting loop containing the 
read-out junction and the chain, enables the control 
of the bias phase 7 = $c — 8 over the chain. 5 is the 
phase difference on the read-out junction. 

We have measured the switching current (see the 
Methods section) of the entire Josephson junction cir- 
cuit containing both the chain and the read-out junc- 
tion. We start by presenting the gate-voltage depen- 
dence of the switching current (Fig. [21 for small vari- 
ations of the gate voltage (so that |Vg| <S 2e/C£) 
and for two different ratios of Ej/Ec ■ As the cen- 
tral coupling Cf is ~ 10 times larger than any other 
coupling, the gate voltage only polarizes the middle 
island. The values of the island charges result from 
the combined effect of the gate voltage and off-set 
charges. In the particular case of our circuit, the latter 
vary randomly within a time scale of ~ 5 min in aver- 
age, enabling the measurement of a single charge con- 
figuration during a gate voltage scan that takes ~ 3 
minutes. Thus, by repeating the same voltage sweep, 
we measure different charge configurations. The re- 
sults presented in Fig. [2] were post-selected from a 
large set of data, of about 3000 runs by choosing the 
largest observed switching current oscillations. The 
largest oscillations displayed in Fig. [2] correspond to 
a particularly simple case in which all charges, ex- 
cept the one on the middle island are close to zero: 

k' = (0± 0.1,0 ±0.1,^,0 ±0.1,0 ±0.1). Notice 
that the ab-initio probability to produce this charge 
configuration is 0.2 4 — 0.0016 which translates into 6 
configurations out of 3000, so its observation supports 
the assumption of random charge distribution. 

The charge frustration on the middle island intro- 
duces a geometrical phase shift exp [i2irq3] between 
the three CQPS occurring on the junctions at the left 
side of the middle island and the three CQPS on the 
junctions at the right of the middle island. This phase 
shift between the different CQPS is graphically rep- 
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Figure 2: Phase-slip interferences controlled by the po- 
larization charge on the middle island of a 6 Josephson 
junction chain, corresponding to the charge configurations 

n' = (0± 0.1,0 ±0.1, ^§^-,0±0.1,0±0.1). (a) Schematic 
representations of the total CQPS amplitude v* (in red), 
obtained by summing up the 6 phase-slip probability am- 
plitudes Vi (in black), represented as vectors in the com- 
plex plane, v* is presented for four different charge con- 
figurations, (b) (c) The black diamonds represent the 
measured variation of the switching current as a function 
of the induced charge on the middle island, in the case 
of Ej/E c = 3 (b) and E.,/E c = 2 (c). The y-axis is 
reported in units of the critical current of one junction, 
io. The red lines represent the theoretical calculations in 
the independent phase slip approximation (MLG model) 
of CQPS interference. The chain was phase biased at a 
constant phase: <!>c = 0.9 n. The working points for the 
measurements presented in (b) and (c) are presented in 
the panels (d) and (e) at the right of each curve. All 
curves are shifted so that zero of the y axis corresponds 
to the switching current of the zero charge configuration 
(0,0,0,0,0). The blue lines represent the corresponding 
calculations using the diagonalization of the Hamiltonian 
dsl . There is a reasonable agreement between the semi- 
classical MLG model of CQPS interference, numerical cal- 
culation and data. 



resented for several charge configurations in Fig. [2^. 
In this regime the total CQPS amplitude (eq. [I]) as a 
function of gate voltage becomes: 



3ua/2 + 2 cos {irVgCl/e). 



(4) 



The phase-slip amplitude is expected to vanish 
completely v* = for the charge configuration 
(0, 0, 0.5, 0, 0) while the maximum value v* = 6v is ob- 
tained for the charge configuration (0,0,0,0,0). The 
red line of Fig. [2] shows the corresponding theoretical 
calculation using the CQPS model (|3]),([4|. Around 
the charge configuration (0,0,0.5,0,0) we expect a 
complete suppression of the total phase-slip amplitude 
v* (see Fig. [2^i), hence an increase of the supercur- 
rent through the chain. For Ej/Ec — 3, the change 
in the measured switching current due to the full sup- 
pression of CQPS is ^lnA, which represents ~ 2% of 
the critical current of one SQUID, io- Increasing the 
CQPS-amplitude by decreasing the ratio Ej/Eq to a 
value of 2 , the oscillation amplitude of the switching 
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current increases to ~ 2nA which represents ~ 5% of 
i (see Fig. (2ji). 

We now turn to the discussion of more complex in- 
terferences of CQPS realized by increasing the gate 
voltage (\Vg\ > 2e/C|) that leads to the polarization 
of the islands next to the central island. In Fig. [3J: we 
show two interference patterns that were post-selected 
from a total of 200 curves. Again, the selection 
criteria was the maximum observed switching cur- 
rent amplitude. In the following we show that these 
measured patterns can be understood by considering 
charge configurations in which only qi,q§ ~ 0: k" — 

(0±0.1,$ + -§^,93 + ~^r~i <W + -fr 1 ^ ±0.1). 
The corresponding ab-initio probability is 0.2 2 = 4% 
which translates into ~ 8 curves out of the measured 
200. 

The charge frustration on the middle, the second 
and the fourth islands introduces geometrical phase 
shifts between the CQPS of the second, third and forth 
junctions. Fig. [3^ show the corresponding CQPS am- 
plitudes as vectors in the complex plane for several 
charge configurations. The resulting switching cur- 
rent oscillations SIsw, presented in Fig. [3J), show a 
complex pattern, composed of a fast harmonic arising 
from the strong Cf coupling and a slower evolving en- 
velope due to the weaker Cf and Cf couplings (see 
Table |l| . In Fig. [3J: we show the measured interfer- 
ence patterns for two different phase biases &c of the 
Josephson junction chain; these biases were chosen 
close to $(7 — t in order to maximize the response 
of the chain. For the top curves in Fig. [3J3 and c 
we polarized the chain at $c £ t so we expect the 
switching current to increase when the phase slips are 
suppressed. Similarly, for the bottom curves, where 
<!>C > 7r, we expect the switching current to decrease 
when the chain becomes classical. The exact shape 
of the oscillations envelope depends on the configura- 
tion of the offset charges q®. For the two calculated 
curves we have chosen the offset charges configura- 
tions , q^ that give the best fit the experimental 
data. The exact values of the fit parameters are shown 
111 table. |U 
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Table I: Fit parameters for the calculated QPS interference 
patterns presented in Fig. [3b 
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Figure 3: CQPS interferences induced by the polariza- 
tion charge on the middle and the first two lateral is- 
lands (0, f/2, <?3, 94, 0) of the 6 Josephson junction chain 
at Ej/Ec = 3. (a) Schematic representations of the 6 
phase-slip probability amplitudes Vi (in black) and the to- 
tal CQPS amplitude v* (in red) as vectors in the com- 
plex plane, for several particularly chosen charge config- 
urations, (b) The calculated switching current oscilla- 
tion, SIsw, induced by the polarization charges for a large 
sweep of the gate voltage Vg, at two different phase biases. 
The polarization charge on the central island 53 is shown 
on the lower x-axis and the charge on one of the lateral 
islands 54 is shown on the higher x-axis. (c) The measured 
SIsw over a large sweep of Vg at the same phase-biases 
$c as in (b). The value of $c for each curve is shown on 
the right side of the figure. 



The qualitative behavior of the interference pat- 
tern agrees perfectly with the theoretical expectations 
based on the simple picture of addition of complex am- 
plitudes. Using the fitted values of q[ ^ we evaluate 
the gate voltages corresponding to the special charge 
configurations shown in Fig. [3^ and indicated them 
on the measured and calculated curves. For instance, 
point (i) corresponds to the configuration where all 
charges are zeros, the same configuration that was re- 
alized in measurements displayed in Fig. [2^ at q$ = 0. 
Notice that in the measurements presented in Fig. [3] 
its position is shifted along the x axis by the offset 
charges on the islands 2 and 4. Because for arbitrary 



values of qi the total CQPS amplitude is less than Nv, 
one expects that maximal oscillations as a function of 
(73 occur around point (i) (see Fig. |3^,). Indeed, next 
to this point, the switching current SIsw oscillations 
have the largest amplitude. The chain goes from the 
perfectly coherent phase-slips regime at (73 = 0, where 
the switching current is minimum (the zero level in 
Fig. [2]d and c), to the maximally dephased configu- 
ration at <?3 = 0.5 when the phase slips are canceled, 
the chain is almost classical and the critical current 
is enhanced. Point (ii) corresponds to opposite limit 
in which the oscillations are strongly suppressed due 
to interference induced by non-zero charges at q^ , qz . 
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Working point (Hi) shows the situation when the to- 
tal CQPS amplitude is suppressed at — and it 
never reaches the maximum Nv for any value of 93 
(see Fig. [3^) . In this case we expect small oscillations 
of SIsw that reach the maximum supercurrent for the 
classical chain. The case of (iv) shows that it is not 
necessary to have the CQPS perfectly aligned as in 
(i) in order to have a large amplitude of SIsw oscilla- 
tions. As expected, the SIsw oscillations around (iv) 
are comparable in amplitude to the ones around point 
(i). 

The Aharonov-Casher interference of phase slips is 
expected to be a 2e periodic effect. In general, ran- 
dom le quasi-particle poisoning is a severe problem 
for the observation of the interference effect as the le 
contamination reduces the accessible charge interval 
from [0, 2e] to [0, le]. As a consequence, the ampli- 
tude of the oscillations in the interference pattern is 
significantly reduced by le quasi-particle poisoning[8 . 
In our case we observe the quasi-particle poisoning by 
the appearance of a le periodicity in the SIsw vs - Vg 
oscillations for temperatures above T — 300 mK. At 
base temperature T — 50 mK we scan the full charge 
space interval [0, 2e] enabling the observation of com- 
plete destructive CQPS interference (see Fig. [2| . Fur- 
thermore, spectroscopy measurements of propagative 
modes in long SQUID chains [TS] independently con- 
firms the value of the gate capacitance Cg ~ 400pF, 
corresponding to a 2e periodic SIsw V8. Vq curve at 
low temperature. 

As can be seen from the periodic dependence of the 
measured switching current, the island charge config- 
uration does not change significantly during the mea- 
surement. Although each measurement point implies 
10 4 repeated switchings into the dissipative state of 
the junctions, where large numbers of quasi-particles 
are excited, after the circuit relaxes back to the dis- 
sipationless state, the charge configuration is stable 
enough in order to enable the measurement of the in- 
terference pattern. We have directly measured the fre- 
quency of random charge jumps by repeating the same 
measurement several times and we observe a typical 
time of T qp ~ 5 minutes between changes in the island 
charge configuration. This time interval is sufficient 
in our case as it enables the measurement of several 
hundreds of experimental data points. 

It is well established that slow drift of charge 
induced by fluctuating TLS leads to Sq <C e at 
time scales of minutes|20 . A significant charge 
drift at this and smaller time scales is attributed to 
non-equilibrium quasi-particles jumping between the 
is lands [2"T|; the equilibration of these quasi-particles 
is made difficult by their localization in subgap 
states|22_. Surprisingly, the time scales of these 
jumps might be dramatically different even in sim- 
ilar devices. Charge fluctuation times, similar to 
the one observed here, have been reported previ- 
ously in small highly resistive Josephson islands[23 , 
in small charge-phase qubits[17| and in Cooper- 
pair transistors |25l 125] . However much shorter times, 
T qp < Is, were reported for the fluxonium circuit [27_ 
and even shorter times, r qp ~ 1/is , were reported for 
larger devices such as the transmon|28|. 



Reasonably good agreement between our measure- 
ments and a model of independent phase slips might 
be surprising given the modest values of Ej/Ec = 
2 — 3 of our Josephson chains, because this approxima- 
tion is expected to be correct only for Ej/Ec 3> 1 [ 1 5J . 
Josephson junction circuits with Ej/Ec ~ 1 are typ- 
ically analyzed using the charge basis description and 
only a couple of charging states are needed for accu- 
rate results. At larger Ej/Eq the number of required 
states grows, which makes the problem of numerical 
diagonalization difficult, especially for large systems. 
From a practical point of view it is important to com- 
pare both approaches and their validity as a function 
of Ej/Ec, because the calculation is many orders of 
magnitude faster in the phase slip approximation, in 
particular for systems containing a large number of 
iunctions|29|. 

We now discuss the details and the validity of the 
diagonalization of the full Hamiltonian of the chain: 

H = 4 4 Ei,i [C~\ (ft - ft) (ft - Qi) + (5) 
+ Ei=l E J I 1 ~ cos (fi -<Pi-l)] 

where Qi is the charge (in units of 2e) on the i-th 
island and <fi is the superconducting phase on the 
island. C^ 1 is the matrix of inverse capacitance of 
the chain. The first sum in the expression (|5j is the 
charging energy for the islands of the chain and the 
second sum represents the Josephson couplings for 
all junctions in the chain. As the total phase differ- 
ence 7 across the chain is fixed, we have (po = and 
^6=7- The charges Qi are multiples of the elemen- 
tary charge of a Cooper pair. As Qi and <p^ are conju- 
gate variables, the chain's wavefunction ip(<pi, 925) 
is 2tt periodic in tp im From the Hamiltonian (|5| one 
can see that its energy spectrum is a periodic func- 
tion of the polarization charges q^. indeed, any mod- 
ification of the polarization charge by the charge of a 
Cooper pair qi — > qi + 1 can be absorbed by the uni- 
tary transformation exp(i^j) which changes Qi into 
Qi + 1, while leaving ipi invariant. Therefore the su- 
percurrent through the chain / (ipi, <fe) remains un- 
changed when the polarization charges change by a 
multiple of the charge of a Cooper pair. 

The blue curves in Fig. and c show the calcu- 
lated switching current oscillations, from the numer- 
ical diagonalization of the total Hamiltonian (|5| (for 
details see the Methods section), as a function of Vq 
for the charge configurations k! . These calculations 
agree reasonably well with the switching current cal- 
culated using the semi-classical CQPS approximation 
(the MLG model) . The modulation of the critical cur- 
rent expected theoretically is somewhat larger than 
the data. This is due to the fact that in the experi- 
ment the random charges qi are not exactly zeros; this 
results in the interference that decreases the observed 
amplitude of the switching current modulation as we 
discussed above for the charge configurations k". 

In conclusion, we have presented a quantitative 
study of the Aharonov-Casher effect exhibited by 
fluxon motion in a multi-junction circuit. We com- 
pared the data with the expectations based on the 
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diagonalization of the full Hamiltonian of the chain in 
the charge basis. Our results show that the ground 
state properties of a short Josephson junction chain 
can be fully understood in terms of phase slip dynam- 
ics even in a parameter range that has been tradition- 
ally described in terms of charge dynamics. The mea- 
surements also show that the polarization charges on 
the islands of the chain can be controlled with suf- 
ficient precision and they are stable enough to en- 
able the observation of the chain's collective behav- 
ior at the time scale of minutes. We believe that 
our results will provide a starting point to reconsider 
the physics of large Josephson junction arrays, long 
Josephson chains and their possible applications to 
the frequency-to-current conversion device or a topo- 
logically protected qubit. 
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Supplementary information: The hopping term in 
the Matveev-Larkin-Glazman theory of quantum 
fluctuations 

Here we present the detailed derivation of the hop- 
ping term v* of the MLG model in the charge frus- 
trated chain. Similar calculations have been per- 
formed for the Josephson chain [15] and for slightly 
different Josephson circuits [301 |3~T] . To calculate the 
hopping term we need to find the classical trajecto- 
ries connecting states before and after one phase-slip 
event. There are N such trajectories, each of them 
corresponding to the phase slip occurring on a par- 
ticular junction in the chain. In a semi classical ap- 
proximation, the contribution of the phase slip in the 
junction i to the hopping term is governed by the 
imaginary-time action Si on the corresponding tra- 
jectory: 

Vi = Ae~ s * (6) 

The prefactor A accounts for the contribution of the 
non-classical trajectories close to the classical one that 
defines Si. 

In order to calculate the actions Si, we need to de- 
rive the complete Lagrangian for the Josephson chain. 
The electrostatic effects in the Josephson chain are de- 
scribed by the following Hamiltonian: 



H ° = |E I ' 1 ] ij(Qi ~ 9i) (Qj - 9i) 



(7) 



Methods 



Switching current measurements The switching 
current was determined from the switching probability 
at 50%. We apply ~ 10 4 bias-current pulses of am- 
plitude Ibias an d measure the switching probability as 
the ratio between the number of switching events and 
the total number of pulses. The measured switching 
current corresponds to the escape process out of the 
total potential energy containing the contributions of 
the read-out junction and the chain. We can calcu- 
late this escape process and therefore deduce the ef- 
fect of quantum phase-slips on the ground state of the 
chain [TS]. From the escape rate, we can deduce the 
escape probability Psw as a function of the current 
bias Ibias and infer the theoretical switching current 
(at Psw — 50%) for each biasing point ($c, Vq)- 

Numerical diagonalization of the exact chain Hamil- 
tonian The Hamiltonian (|5| gives the exact descrip- 
tion of the Josephson circuit (provided that all en- 
ergy scales remain small compared to Cooper gap A 
and no other degrees of freedom are involved. For 
the purposes of numerical diagonalization one has to 
limit the number of charging states on each island. 
This approximation can be easily tested by compar- 
ing the results of diagonalization for different number 
of allowed charging states. For the problem here with 
Ej/Ec = 2 — 3 it is sufficient to keep 7 charging states 
to get the results with 10 -2 accuracy. 



The polarization charges qi 



CfV a 



axe controlled 



by the gate voltage. We would like to mention that 
in our experimental setup we have added screening 
lines to the central gate, in order to obtain a coupling 
to the central island at least 10 times larger than the 
couplings to the rest of the chain: Cf ~ 10 * Cf, Cf — 
50*Cf,Cf. 

Since the charges Qi and the phases of the islands 
tfii are canonical conjugate variables, the equation of 
motion for the phase reads: 



*-^-E[c-V<>i-*> 



J J 

(8) 

Using eq. ^ we can rewrite the charging Hamilto- 
nian ([7]) in the phase notation: 



ijfifj 



(9) 



The charge part of the Lagrangian for the Josephson 
junction chain reads: 



(10) 



Following formula ( 10 ) and using the expressions (IS]) 
and ([9]) we get for the charge Lagrangian the following 
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expression: 

Cc = E CiiVWj + E ~ \ E CiWWj 



4Eg 



E 



q%<P% 



(11) 



The capacitance matrix Cy contains the values of 
all coupling between the islands. However, in real- 
ity, due to the geometry of the sample, the capaci- 
tance between first neighbors is orders of magnitude 
larger then the stray capacitance between second or- 
der neighbors. This means that we can safely work 
within the so called nearest neighbor capacitance ap- 
proximation, and the matrix dj only gets non zero 
contributions for the elements closest to the main di- 
agonal: 



/ 2C -C 
-C 2C -C 
-C 2C 

V 



... ^ 

... 1 

... 

... -C 

-C 2C J 



(12) 



Where C is the capacitance of one junction in the 
chain. 



Using the approximation ( 12 1 the expression of the 
charge Lagrangian is simplified and it reads: 



Cr 



<pi-i) + 



E 



(13) 



Introducing the phase differences on the junctions 
0% = fi+i — <Pi and including the Josephson energy, 
we derive the complete Lagrangian of the chain: 



E 



16E C 



— Ej cos 



We can see that the Lagrangian ( 14 1 has two com- 



i-l 

£* 

(14) 



ponents which have very different physical conse- 
quences. The first sum that we call Co is inde- 
pendent on the frustration charges q^. It gives a 
contribution to the real part of the phase-slip am- 
plitude Vi, that is given by the Bloch band width 



0.2F, 



e V E c . For identical junc- 
tions in the chain, v is independent on the path chosen 
by the phase slip. The second sum of the Lagrangian 
( fl4| , which we call S£, has the form of a total time 
derivative. Hence, this term does not change the clas- 
sical equations of motion and the real part of the clas- 
sical action on a single trajectory. However, 8C gives 
the tunneling amplitude along each path its own phase 
factor. This phenomenon is mathematically equiva- 
lent to the AB effect for the phase variable ifi which is 



2n periodic. Changing pi amounts to changing the pe- 
riodic boundary conditions ip(cpi + 2tt) = e l27Tpi ip((pi) 
for the phase, in analogy to the motion of a charged 
particle on a circle, threaded by a flux tube. 

When a phase-slip occurs on junction i, the other 
phase differences 6j are changed by: 



A8j = J?L 
N 



2ttS, 



(15) 



Thus, the contribution to the phase-slip action from 
the j-th junction in the presence of charge frustration 
reads: 



SSj = —i / SCdt = — i^^pfeA^fc = —2nipj 



2ni 



N ^ Pk 

k 

(16) 

Since the last term in the expression above does not 
depend on k, it only adds an overall phase term for 
all phase slip trajectories, thus has no physical effect 
on the interference pattern and it can be dropped. 
Replacing this result in the formula (|6J, we get the 
mathematical expression for the charge frustration de- 
phasing factor in the phase slip probability amplitude 
of the j-th junction: 



(17) 



So the phase slip probability amplitude on the j-th 
junction vj reads: 



v exp 



t2nJ2 



Qk 



fc=i 



(18) 



In other words, the absolute value of the probability 
amplitude for the QPS is the same as in the absence of 
charge frustration, but the geometric phase difference 
between the QPS is proportional to the total charge 
on the islands between the junctions. Finally, the full 
hopping term between the states |m) and \m + 1) in 
the presence of charge frustration is the sum of phase 
slip amplitudes m in all six junctions: 



E'' 

i=l 



(19) 



At zero gate voltage, the expression (19 1 reduces 
to v* = Nv that was used in the previous section to 
solve the tight binding MLG Hamiltonian and calcu- 
late the expected switching current. Non-zero gate 
voltage directly affects the interference of QPS by 
changing the geometrical Aharonov-Casher phase dif- 
ference between phase slips in different junctions and 
thus provides a direct test for the quantum nature of 
the chain's ground state. 
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